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ABSTRACT 


The  principle  of  dynamic  equipollency  is  defined 
and  applied  to  two  areas  of  mechanical  admittance/ 
impedance  testing:  generation  of  experimental 

moment -rotational  velocity  admittances  and  im- 
pedances, and  many-DOF  structures. 
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APPLICATIONS  OF  THE  PRINCIPLE  OF  DYNAMIC 
EQUIPOLLENCY  TO  MECHANICAL  ADMITTANCE  TESTING 


INTRODUCTION 

The  principle  of  dynamic  equipollency  (DE)  is  the  relationship  between 
two  sets  of  forces  applied  at  two  separate  subsets  of  a structure's  degrees  of 

freedom  (DOF).  This  report  illustrates  and  elucidates  applications  of  the  DE 
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principle,  as  proposed  in  Lincoln  Laboratory  Technical  Note  1974-17,  to  two 
weak  areas  of  mechanical  admittance/impedance  techniques:  (a)  many  DOF 

structures,  and  (b)  experimental  generation  of  moment-rotational  velocity 
impedances. 

The  DE  principle  is  applied  to  a structure  with  25  DOF  to  determine  the 
structure's  first  two  mode  shapes  and  resonant  frequencies,  greatly  reducing 
the  number  of  measured  mechanical  admittances  required  to  characterize  a 
many -DOF  structure. 

Also,  moment -rotational  velocity  impedances,  determined  from  measured 
admittances  via  the  DE  principle,  make  possible  the  assembly  of  substructures 
into  a total  structure  ensuring  the  compatibility  of  rotational  as  well  as  lateral 
displacements  at  interfaces.  Two  beam  structures  are  assembled  mathemat- 
ically using  experimentally  determined  impedances. 

It  is  important  that  the  reader  have  a basic  understanding  of  mechanical 
admittance/impedance  concepts.  It  is  also  desirable  for  the  reader  to  be  fa- 
miliar with  the  inverse  eigenvalue  method  of  mechanical  admittance  analysis, 
a method  of  solving  for  a structure's  mass  and  stiffness  matrices  from  me- 

3 

chanical  admittance  measurements. 
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CHAPTER  I 

PRINCIPLE  OF  DYNAMIC  EQUIPOLLENCY 


When  the  number  of  mode  shapes  of  a large  complex  structure  is  less  than 
the  number  of  locations  on  the  structure  (characterized  as  DOF)  where  dynamic 
behavior  is  important,  certain  problems  are  encountered  with  the  structure's 
admittance  matrix.  The  order  of  the  admittance  matrix  is  equal  to  the  number 
of  DOF  where  admittance  measurements  are  made  and  the  rank  of  the  admit- 
tance matrix  is  equal  to  the  number  of  mode  shapes  (rigid  body  and  elastic). 

In  this  case,  the  order  of  the  admittance  matrix  is  greater  than  the  rank  of 
the  admittance  matrix.  Consequently,  the  admittance  matrix  is  singular  and 
its  inverse,  the  impedance  matrix,  does  not  exist. 

The  columns  (rows)  of  this  singular  admittance  matrix  are  linear  combina- 
tions of  fewer  mode  shapes  than  there  are  columns  (rows)  in  the  matrix.  It 
follows  that  there  must  exist  relationships  among  the  columns  (rows)  of  the 
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admittance  matrix.  These  relationships  are  determined  using  a linear,  modal 
model  of  the  structure. 

Equation  (1)  is  an  expression  of  the  velocities  at  all  of  the  structure's  DOF, 
{v},  resulting  from  the  applied  forces  given  in  the  column  matrix  {f}.  L,  the 
number  of  DOF,  is  greater  than  X,  the  number  of  mode  shapes  of  the  structure: 


X =£-  U'J  » 

Y I 


(1) 


i=1  -wi)+D/a;) 


where 


{$*}  = £ ^ mode  shape  of  the  structure,  { = 1,  2,  3.  . . X 
mi  = * ^ modal  mass 

oo 2 = i th  resonant  frequency  expressed  in  radians  per  second 


D^(o;)  = f ^h  modal  damping  function. 
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The  column  matrices  {0^},  {f}  and  {v}  have  L elements.  Each  column 
matrix  is  partitioned  into  two  subcolumn  matrices  corresponding  to  two  mutu- 
ally exclusive,  collectively  exhaustive  subsets  of  the  structure's  DOF.  The 
first  subset  is  the  DOF  f = 1,  2,  3.  . . and  the  second  subset  is  the  DOF 
t = X + 1,  \ + 2,  X + 3,  . . . L. 


{/}  = 


W = 


(f)  = 


(2) 

(3) 

(4) 


Equation  (5)  is  Eq.  (1)  expressed  using  Eqs.  (2),  (3)  and  (4): 

[i 

f2 


•41  - 2 


l HT  KkzJ  14 


f =1 


i(aj2  - oj  2)  + D (a; ) 


0, 


0- 


= z 


1=1 


00  i j.  4 . ^ | I fl  I 

m^  L^l  \^z\  I O'  I 

i(aj2  - oj  2)  + D (a; ) 


0. 

0. 


V1  I _ V 

V2  1 2 , . 

f =1 


•v  I l 1 I , I I I u 

* mT  L01  ; 02 J If, 


i(co2  - a;  2)  + D^oj) 


*1 

77 

^2 


(5) 


(6) 


(7) 


Equations  (6)  and  (7)  represent  cases  where  forces  are  applied  only  at  the 
first  DOF  subset,  and  only  at  the  second  DOF  subset,  respectively.  If  the 
relationship  in  Eq.  (8)  is  satisfied,  the  dynamic  response  of  the  structure  is 
the  same  for  the  two  cases  shown  in  Eqs.  (6)  and  (7): 


i = 1,  2,  3.  . . X . (8) 
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Equation  (9)  is  Eq.  (8)  simplified: 

[0lJ  V = 1^2 J (f2}  f = 1,  2,  3. . . \ . (9) 

The  set  of  X equations  shown  in  Eq.  (9)  are  consolidated  into  one  matrix  equa- 
tion in  Eq.  (10). 

[01JT  {f4>  = [<^2]T  (f2>  • (10) 

T 

The  matrix  [$~]  is,  in  general,  rectangular  and  has  rows  that  are  the  struc- 

^ T 

ture's  mode  shapes  at  the  second  DOF  subset.  [</>,]  is  a square  matrix  and 

1 q1 

its  rows  are  the  structure's  mode  shapes  at  the  first  DOF  subset.  Since  [<t>^] 
is  square  and  its  rows  are  X,  linearly  independent  mode  shapes,  then 
1 exists: 


{f4}  = U1]T"1  [4>2]T{f2} 


(11) 


Equation  (11)  states  that  given  any  set  of  forces  applied  at  the  second  DOF 
subset,  {f ^ },  one  can  solve  for  a set  of  forces  to  be  applied  at  the  first  DOF 
subset,  {f^},  that  produces  the  same  dynamic  response  of  the  structure.  This 
is  the  DE  principle. 

The  structure's  admittance  matrix  is  partitioned  in  Eq.  (12)  according  to 
the  previously  mentioned  DOF  subsets: 


Xi. 

[Yll 

Y1  2 

V2 

Y21 

Y22 

(12) 


|V1 

r Yi  i ! Yi  2 1 

|fl 

>V2 

1 

.Y21:Y22. 

1 0 

1 Vl  1 

YH  i Y12  1 

1 0 

1 v2 

1 2 

Y2! ! Y22  J 

' "f2 

(13) 


(14) 


It  has  already  been  shown  that  the  two  sets  of  applied  forces  {f^}  and  {f^|  can 
be  determined  such  that  they  produce  the  same  dynamic  response  of  the  struc- 
ture. Equations  (13)  and  (14)  can  now  be  used  to  determine  relationships  among 
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the  columns  (rows)  of  the  structure's  admittance  matrix.  The  columns  (rows) 
are  expressed  collectively  as  submatrices: 


y : y 
*1 1 J_ 1 1 2 

i.T 

= I-1- 

_ |VH 

Y ' Y 
_ llj  _ 12 

Y"  ' Y 

21  ! 22  J 

1 o 

V2 

1-^ 

J 

<! 

! 

DJ 

1 

V ' V 

X 21 ! *22  J 

[Y11){f1}=  IY12J  {£2} 

[Y21]{f1}=  (Y22]{f2)  . 

Equation  (19),  derived  from  Eq.  (17),  is  similar  to  Eq.  (11): 

<fl>  = tY11l_1  [Y12l  {£.,} 

{q>  = [01JT‘1  (4>2IT  (r2) 


(15) 

(16) 

(17) 

(18) 

(19) 

(11) 


[Y12]=[Y21]T  • (20) 

Equations  (11)  and  (19)  and  the  symmetry  of  the  admittance  matrix,  Eq.  (20), 
relate  measured  admittances  to  the  structure's  mode  shapes,  Eq.  (21): 

[01] T“1  [<PZ]T  = [Y^f1  t Y21)T  • (21) 

“I 

The  matrix  [Y^J  can  be  calculated  from  [ ^ J except  when  [ ^ J is  meas- 

ured at  a resonant  frequency  of  the  structure. 

-1  T 

Equation  (21)  shows  that  the  matrix  product  [Y^j  [Y^]  is  invariant 
with  respect  to  the  frequency  of  excitation.  It  is  important  to  remember  that 
one  of  the  assumptions  leading  to  Eq.  (21)  is  that  the  structure  has  only  a finite 
number,  X,  of  mode  shapes.  In  general,  the  elements  of  [ Y^ ^ ] and  [ Y^ ^ ] are 
functions  of  the  frequency  of  excitation.  Equation  (21)  yields  Eq.  (22)  and  is 
used  in  Chap.  II  to  solve  for  the  mode  shapes  of  a 25 -DOF  beam  structure: 

[02] T  = [01]T  [Y11]“1  [Y21]T  . (22) 
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The  DE  principle  can  be  used  to  determine  the  relationship  among  all  the 
submatrices  of  the  structure's  admittance  matrix.  Substituting  Eq.  (19)  into 
Eq.  (18)  gives  Eq.  (23): 


^ Y2iJ  ^Yll^  ^Y12^  lf2'  ^ Y zn}  ‘^2 1 


(23) 


Equation  (23)  yields  Eq.  (24)  which  is  the  relationship  among  the  submatrices 
of  the  admittance  matrix: 


[y2iHYn 


-l 


[y12I  =[y22] 


(24) 


Equation  (25)  is  derived  using  the  symmetry  of  the  structure's  admittance 
matrix,  Eq.  (20),  and  the  relationship  in  Eq.  (24): 


Y ) • f Y 1 • 1 Y 1 Y 

1 21  1 11J  1 21J 


lY22l 


(L-X)XX  XXX  XX(L-X)  (L-X)X(L-X) 


(25) 


The  number  of  elements  in  each  matrix  is  given  under  each  matrix. 

As  before,  X is  the  number  of  modes  of  the  structure  and  L is  the  number 
of  DOF  for  the  structure  where  L > X.  If  L is  much  greater  than  X,  it  can  be 
seen  from  Eq.  (25)  that  it  is  much  more  efficient  to  measure  the  elements  of 
[ Yf  i J and  [Y,,^]  (XXL  total  elements)  and  calculate  [ Y^J  using  Eq.  (25)  rather 
than  measure  the  (L  — X)^  total  elements  of  [Y^]. 

Equation  (25)  can  also  be  used  to  calculate  moment-rotational  velocity 
admittances.  If  the  second  DOF  subset  is  rotational  DOF  and  [Y^]  is  meas- 
ured using  rotational  sensors  at  the  second  DOF  subset,  [Y^J  calculated  using 
Eq.  (25)  is  comprised  of  moment-rotational  velocity  admittances: 


iy21i 


lyU> 


-1 


[ Y 1 

1 I21 1 


^22^ 


/i  I v “1  /lb-sec,  X-1  ,,,  . x “1 

(lb-sec)  • ( — : ) • (lb-sec)  = (lb-in. -sec) 


in. 


(25) 

(26) 


The  physical  dimensions  of  each  matrix  in  Eq.  (25)  are  given  directly  under  it 
in  Eq.  (26).  Equation  (26)  shows  that  [Y^J  has  the  physical  dimensions  of 
moment-rotational  velocity  admittances.  The  experimental  determination  and 
application  of  moment-rotational  velocity  admittances  are  the  subjects  of 
Chaps.  Ill  and  IV. 
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CHAPTER  II 

DETERMINATION  OF  MODE  SHAPES  FOR  A 25 -DOF  STRUCTURE 


A 6 -ft,  clamped -clamped  beam  structure  is  characterized  as  having 
25  DOF,  L = 25  (Fig.  1).  The  aluminum  channel  beam  is  clamped  at  both  ends 
to  a steel  support.  The  structure  is  assumed  to  have  two  mode  shapes,  X = 2, 
that  dominate  the  structure's  low  frequency  (less  than  150  Hz)  dynamic  behavior. 
The  first  DOF  subset  is  DOFs  10  and  16.  The  second  DOF  subset  consists  of 
all  the  remaining  DOFs. 


18-7-2717 


MASS  PER  UNIT  LENGTH  3.74  X lo"4sluq/in. 

LENGTH  72.0  in. 

Fig.  1.  Six-ft,  clamped-clamped,  25 -DOF  beam  structure. 
All  DOF  are  2 -in.  apart. 


The  DOFs  on  the  structure  are  chosen  such  that  no  axial  or  tortional  motion 
of  the  beam  is  generated  while  the  admittances  are  measured.  The  Wilcoxon 
Research  F3  shaker  and  Z602  impedance  head  assembly  generate  and  measure 

5 

the  applied  force  and  resultant  motion  of  the  structure.  The  admittances  are 
generated  by  a Spectral  Dynamics  Corporation  SD-1002B-33  automatic  me- 
chanical impedance  system  (Fig.  2)/*  A piezoelectric  accelerometer  measures 
motion  of  the  structure  at  DOFs  other  than  the  one  where  the  force  is  applied. 
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Fig.  2.  SD1002B-33  automatic  mechanical  impedance  system  (left) 

and  power  and  charge  amplifiers  (right). 
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The  shaker/ impedance  head  is  placed  at  DOF  16  and  the  mutual  admittance 
and  all  cross  admittances  are  measured.  Since  DOFs  10  and  16  are  equidistant 
from  the  center  of  the  beam  and  the  beam  is  symmetric  about  its  center,  the 
motion  of  the  beam  with  the  shaker/ impedance  head  attached  at  DOF  10  is  a 
mirror  image  of  the  motion  of  the  structure  with  the  shaker/impedance  head 
attached  at  DOF  16.  The  admittances  measured  with  the  shaker/impedance 
head  at  DOF  16  are  converted  to  admittances  of  the  structure  with  the  shaker/ 
impedance  head  attached  at  DOF  10  using  the  scheme  shown  in  Eq.  (27): 

YS  16  = Y26-S  10  • (27* 

Values  chosen  from  the  mutual  and  cross  admittances  measured  for  the 
first  DOF  subset,  DOFs  10  and  16,  are  used  to  calculate  mass  and  stiffness 
matrices  of  the  structure.  The  mass  and  stiffness  matrices  are  calculated 
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using  the  inverse  eigenvalue  method: 


[M] 


' 8.75 

.-3.29 


-3.29 

8.75. 


slug 


(28) 


[K] 


4.89  -4.28 

.-4.28  4.89 


102  lb/in. 


(29) 


The  mode  shapes  for  DOFs  10  and  16  and  the  resonant  frequencies  of  the  struc- 
ture are  calculated  from  the  mass  and  stiffness  matrices  using  the  iterative 

g 

method.  The  resonant  frequencies  calculated  from  the  experimentally  deter- 
mined mass  and  stiffness  matrices  are  in  good  agreement  with  the  resonant 
peaks  observed  in  the  measured  admittances.  This  lends  credence  to  the 
accuracy  of  the  experimentally  determined  mass  and  stiffness  matrices: 


{ <P 


1.0 

1.0 


1 

53.2  Hz 

1 2 7T 

U1 

= 139  Hz 

2 2? r 


(30) 


(31) 
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Equation  (22)  is  used  to  calculate  the  components  of  the  structure's  mode 
shapes  at  the  second  DOF  subset. 

[4>2]T  = [4>1JT  [Y^r1  [Y21]T  . (22) 

T T 

[4>^]  is  constructed  from  Eqs.  (30)  and  (31).  [Y^J  anc*  [Y^J  are  constructed 

from  the  admittances  measured  at  70  Hz.  The  subscripts  of  the  measured  ad- 
mittances show  the  position  of  the  measured  admittances  in  matrices  [Y^]  and 
[Y.>1]r^  in  Eqs.  (33)  and  (34): 


[< P 4J 


T 


[ Yi  i J 


1.0  1.0' 
1.0  -1.0. 


'4.61 

5.29' 

-1 

10, 10 

10,  16 

.5.29 

4.61  . 

• (-10  ) i 

.16,  10 

16,  1 6 . 

(32) 


(33) 


[y21] 


1.5  1.8  2.3  2.6  3.1  3.5  3.8  4.1  4.3  4.7  5.1  5.3  5.3  5.3 
2.2  2.6  3.2  3.6  4.1  4.4  4.8  5.1  5.3  5.3  5.3  5.3  5.1  4.7 


5.3  5.1  4.8  4.4  4.1  3.6  3.2  2.6  2.2' 

4.3  4.1  3.8  3.5  3.1  2.6  2.3  1.8  1.5. 


(-1 0_1 ) i 


'1,10  2,  10  3, 10  4,  10  5,  10  6,  10  7,  10  8,  10  9,  10  11,10 
.1,16  2,  1 6 3,16  4,16  5,16  6,16  7,16  8,  16  9,16  11,16 


12,  10  13,10  14,  10  15,  10  17,  10  1 8,  10  19,  10  20,  10 
12,  16  13,16  14,  16  15,  16  17,  16  18,  16  19,  16  20,  16 

21,10  22,10  23,10  24,10  25,  10' 

21, 16  22,  16  23,  16  24,  16  25,  16  . 


(34) 
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— ^ rp  rp  rp 

[Y^]  is  calculated  and  then  [(/^J  i-s  calculated  from  , 1^21^  an<^ 

[Y^J  , Eq.  (22).  The  calculated  mode  shapes  for  the  clamped-clamped  beam 
structure,  normalized  by  the  largest  component  of  the  mode,  are  compared 
to  the  modes  of  a clamped-clamped  continuum  beam  (Figs.  3,  4 and  Table  1). 

(See  the  Appendix  for  methods  used  to  calculate  the  mode  shapes  of  the  clamped- 
clamped  continuum  beam.) 

The  values  calculated  for  the  second-mode  shapes,  {$  },  are  less  accurate 
than  those  values  calculated  for  the  first -mode  shape,  {(/>*}.  The  frequency 
of  excitation  for  the  values  of  [ ^ ] and  [Y^]  used  to  calculate  [0^]  is  ^0  Hz. 
This  frequency  is  relatively  close  to  the  first  resonant  frequency  of  the  struc- 
ture, 50  Hz,  and  relatively  far  from  the  second  resonant  frequency,  138  Hz. 
Consequently,  the  participation  of  the  first  mode  is  much  greater  than  partici- 
pation of  the  second  in  the  admittance  matrices  [Y^]  and 

[ Yf  i j and  [Y^]  are  measured  at  a frequency  of  excitation,  126  Hz,  rela- 
tively close  to  the  second  resonant  frequency  of  the  structure.  [Y^]  and 
[Y^]^1  measured  at  126  Hz  are  given  in  Eqs.  (35)  and  (36): 


[Ylli  = 


0.93  -3.72 

-3.72  0.93 J 


-1 

10  i 


(35) 


[Y2i] 


1.8  2.5  2.7  2.8  2.9  2.7  2.5  2.3  1.9  0.6 

-3.5  -4.1  -4.9  -5.0  -5.5  -6.0  -5.4  -4.8  -4.6  -3.3 

-0.5  -1.6  -2.5  -3.3  -4.6  -4.8  -5.4  -6.0  -5.5  -5.0 

-2.5  -1.6  -0.5  0.6  1.9  2.3  2.5  2.7  2.9  2.8 


-4.9  -4.1  -3.5 
2.7  2.5  1 . 8 J 


10_1  i 


(36) 


-1  T 

Again  [Y.  is  calculated  and  Eq.  (22)  is  used  to  calculate  [<P^  • The  mode 

shapes  determined  from  the  admittances  in  Eqs.  (35)  and  (36)  are  compared  to 
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MAGNITUDE  NORMALIZED  BY  LARGEST  COMPONENT 
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Fig.  3.  Experimentally  determined  first  mode  shape  for  the  6 -ft 
beam  structure  compared  with  the  first  mode  shape  for  a clamped- 
clamped  continuum  beam. 
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Fig.  4.  Experimentally  determined  second  mode  shape  for  the  6 -ft 
beam  structure  compared  with  the  second  mode  shape  for  a clamped- 
clamped  continuum  beam. 
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TABLE  1 


FIRST  AND  SECOND  MODE  DATA 
FOR  ADMITTANCES  MEASURED  AT  70  Hz 


DOF 

First  Mode 

Second  Mode 

Experiment 

Continuum 

Beam 

Experiment 

Continuum 

Beam 

1 

0.33 

0.30 

0.7 

0.66 

2 

0.42 

0.37 

0.8 

0.77 

3 

0.52 

0.46 

0.9 

0.89 

4 

0.58 

0.54 

1.0 

0.96 

5 

0.68 

0.63 

1.0 

0.99 

6 

0.75 

0.71 

0.9 

1.00 

7 

0.81 

0.78 

1.0 

0.95 

8 

0.87 

0.85 

1.0 

0.87 

9 

0.91 

0.91 

1.0 

0.75 

10 

0.93 

0.95 

0.7 

0.58 

11 

0.94 

0.98 

0.6 

0.40 

12 

0.98 

1.00 

0.2 

0.23 

13 

1.00 

1.00 

0.0 

0.00 

14 

0.98 

1.00 

-0.2 

-0.23 

15 

0.94 

0.98 

-0.6 

-0.40 

1 6 

0.93 

0.95 

-0.7 

-0.58 

17 

0.91 

0.91 

-1.0 

-0.75 

18 

0.87 

0.85 

-1.0 

-0.87 

19 

0.81 

0.78 

-1.0 

-0.95 

20 

0.75 

0.71 

-0.9 

-1.00 

21 

0.68 

0.63 

-1.0 

-0.99 

22 

0.58 

0.54 

-1.0 

-0.96 

23 

0.52 

0.46 

-0.9 

-0.89 

24 

0.42 

0.37 

-0.8 

-0.77 

25 

0.33 

0.30 

-0.7 

— 0.66 
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Fig.  5.  Experimentally  determined  first  mode  shape  for  the  6 -ft 
beam  structure  compared  with  the  first  mode  shape  for  a clamped- 
clamped  continuum  beam. 
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Fig.  6.  Experimentally  determined  second  mode  shape  for  the  6 -ft 
beam  structure  compared  with  the  second  mode  shape  for  a clamped- 
clamped  continuum  beam. 
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the  mode  shapes  of  a clamped-clamped  continuum  beam  (Figs.  5,  6 and  Table  2). 
The  second -mode  shape  is  more  accurate  than  the  first-mode  shape  because 
the  matrices  [Y^j  and  [Y^j  measured  at  126  Hz  are  dominated  by  participation 
of  the  second  mode. 

One  must  ensure  that  participation  of  the  mode  shape  to  be  calculated  using 
Eq.  (22)  is  large  compared  to  the  participation  of  the  other  mode  shapes  in 
[ Yf  i ] and  [Y^].  As  a result,  each  mode  shape  of  the  structure  should  be  cal- 
culated from  values  of  [Y^]  and  [Y^J  measured  at  a frequency  of  excitation 
relatively  close  to  the  mode  shape's  corresponding  resonant  frequency.  The 
more  accurate  first-mode  shape  (Fig.  3)  is  calculated  from  data  measured  at 
70  Hz,  relatively  close  to  the  first  resonant  frequency  observed  in  the  meas- 
ured admittances,  50  Hz.  The  more  accurate  second-mode  shape  (Fig.  6)  is 
calculated  from  data  measured  at  126  Hz,  relatively  close  to  the  second  reso- 
nant frequency  observed  in  the  measured  admittances,  138  Hz.  However, 
choosing  values  of  [Y^j  and  [Y^J  at  frequencies  too  close  to  the  resonant 
frequencies  can  result  in  inaccurate  mode  shapes. 

At  and  near  resonant  frequencies  of  the  structure,  the  mechanical  admit- 
tances are  governed  to  a large  degree  by  the  damping.  If  the  damping  is  signif- 
icantly nonlinear,  the  mechanical  admittances  measured  near  resonant  fre- 
quencies of  the  structure  are  not  of  the  form  shown  in  Eq.  (37): 


If  an  admittance  of  the  structure,  measured  near  a resonant  frequ-ency,  cannot 
be  modeled  by  Eq.  (37),  then  the  response  of  the  structure  cannot  be  modeled 
by  Eq.  (1)  and  the  DE  principle  does  not  hold.  If  at  frequencies  of  excitation 
away  from  resonant  frequencies  the  contribution  of  the  damping  is  relatively 
small,  Eq.  (3  8),  the  dynamic  response  of  the  structure  can  be  modeled  by 
Eq.  (1)  and  the  DE  principle  can  be  applied: 


to 


(37) 


f =1  v t 1 r ' 

Equation  (37)  is  the  expression  for  the  nj^1  admittance  derived  from  Eq.  (1). 


i = 1,  2,  3.  . . A 


(38) 
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TABLE  2 


FIRST  AND  SECOND  MODE  DATA 
FOR  ADMITTANCES  MEASURED  AT  126  Hz 


First  Mode 

Second  Mode 

DOF 

Experiment 

Continuum 

Beam 

Experiment 

Continuum 

Beam 

1 

0.51 

0.30 

0.61 

0.66 

2 

0.48 

0.37 

0.76 

0.77 

3 

0.66 

0.46 

0.87 

0.89 

4 

0.66 

0.54 

0.89 

0.96 

5 

0.78 

0.63 

0.97 

0.99 

6 

1.00 

0.71 

1.00 

1.00 

7 

0.87 

0.78 

0.91 

0.95 

8 

0.75 

0.85 

0.81 

0.87 

9 

0.81 

0.91 

0.75 

0.75 

10 

0.84 

0.95 

0.55 

0.58 

11 

0.81 

0.98 

0.44 

0.40 

12 

0.90 

1.00 

0.23 

0.23 

13 

0.93 

1.00 

0.00 

0.00 

14 

0.90 

1.00 

-0.23 

-0.23 

15 

0.81 

0.98 

-0.44 

-0.40 

1 6 

0.84 

0.95 

-0.55 

-0.58 

17 

0.81 

0.91 

-0.75 

-0.75 

18 

0.75 

0.85 

-0.81 

-0.87 

19 

0.87 

0.78 

-0.91 

-0.95 

20 

1.00 

0.71 

-1.00 

-1.00 

21 

0.78 

0.63 

-0.97 

-0.99 

22 

0.66 

0.54 

-0.89 

-0.96 

23 

0.66 

0.46 

-0.87 

-0.89 

24 

0.48 

0.37 

-0.76 

-0.77 

25 

0.51 

0.30 

-0.61 

-0.66 
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If  the  damping  is  large  and  [Y^]  and  [Y^J  cannot  be  measured  at  fre- 
quencies of  excitation  close  enough  to  the  resonant  frequencies  to  calculate 
accurate  mode  shapes  without  the  damping  significantly  distorting  the  admit- 
tances, there  is  an  alternative. 

The  DE  principle  can  be  applied  with  less  stringent  restraints  than  Eq.  (38) 
by  ignoring  completely  the  real  part  of  the  measured  admittances.  Equation  (39) 
is  Eq.  (37)  rewritten  to  separate  the  real  and  imaginary  parts  of  the  admittance: 


Y . = 
nJ 


■v  CO  •/  2 2.  , _ , 

X ~wf)+D/w)] 


Z 

1=1 


(CO2  -co2)2  + D2(co) 


Im(  Y •)  = 
nj' 


Z 

1=1 


co  t t . , 2 2 

<t>„  <t>-  i(co  . -CO  ) 


n J 


2 2 2 2 
(w  — co  + D (co) 


(39) 


(40) 


Equation  (40)  is  the  imaginary  part  of  the  admittance.  If  the  conditions  in 
Eq.  (41)  are  met  then  Eq.  (42)  can  be  derived  from  Eq.  (40): 


(co2  - co2)2  » D2(co) 


= 1,  2,  3.  . . \ 


(41) 


Im<Ynl)  * Z 


x ASL*V 

n J 


, (co  " — co") 
1=1  w > 


(42) 


The  conditions  imposed  on  the  choice  frequencies  of  excitation  where  the 
admittances  should  be  measured  in  Eq.  (41)  are  less  restrictive  than  the  condi- 
tions imposed  by  Eq.  (38).  For  instance,  with  the  same  damping,  D^(co),  the 
conditions  in  Eq.  (41)  allow  one  to  choose  frequencies  of  excitation,  co,  closer 
to  the  resonant  frequencies  of  the  structure  than  does  Eq.  (3  8).  Also  for  the 
same  frequency  of  excitation,  co,  Eq.  (41)  allows  the  damping,  D^(co),  to  be 
larger  than  does  Eq.  (38). 

Equation  (42)  indicates  that  as  long  as  the  conditions  in  Eq.  (41)  are  met, 
the  damping  no  longer  significantly  affects  the  imaginary  part  of  the  measured 
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admittances.  As  a result,  the  imaginary  parts  of  measured  admittances  ap- 
proximate the  undamped  admittances  of  the  structure: 


Only  the  imaginary  parts  of  the  measured  admittances  are  used  in  calculations 
in  this  report. 

Mechanical  admittances  measured  close  to  resonant  frequencies  of  the 
structure  can  be  inaccurate  for  reasons  other  than  damping.  The  mass  of  the 
shaker,  the  additional  damping  introduced  by  the  shaker  and  the  change  in  stiff- 
ness of  the  structure  due  to  the  attachment  of  the  shaker  have  local  maxima 

9 

of  their  effects  on  measured  admittances  at  resonant  frequencies.  In  addition, 
if  the  structure  is  suspended  to  simulate  a free  structure,  the  impedance  of 

the  suspension  system  has  local  maxima  of  its  effects  on  measured  admittances 

10 

at  resonant  frequencies  of  the  structure. 

-1  T 

By  comparing  the  value  of  [ Y.  . ] [ . J at  70  Hz  vs  the  value  of 

-IT  11  -i  t 

[Y^j  [ Y^i  J at  126  Hz  it  is  apparent  that  [Y^]  [Y^]  is  not  invariant 

with  respect  to  the  frequency  of  excitation.  This  is  because  the  structure  has 
more  than  two  mode  shapes  and  corresponding  resonant  frequencies.  The 
higher  the  frequency  of  excitation  at  which  [Y^j  and  [Y^]  are  measured,  the 
larger  the  participation  of  third  and  higher  mode  shapes  in  the  measured  ad- 
mittances. By  modeling  the  structure  with  more  than  two  mode  shapes, 

-1  T 

[ Y^j  [ i ] should  be  less  sensitive  to  changes  in  the  frequency  of  excitation 
within  the  0—150  Hz  range. 

This  method  of  experimentally  determining  a structure's  mode  shapes  has 
certain  advantages  over  methods  that  attempt  to  excite  and  measure  one  mode 
shape  of  a structure  at  a time.  First,  less  sophisticated  equipment  is  required 
to  measure  an  admittance  than  is  necessary  to  produce  a set  of  applied  forces 
orthogonal  to  all  but  one  mode  shape  of  a structure.  In  addition,  it  is  difficult 

to  produce  a set  of  applied  forces  that  is  orthogonal  to  all  but  one  quasi-known 

1 1 

mode  shape  of  a structure.  Usually  more  than  one  mode  shape  is  contributing 


(43) 
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to  the  dynamic  response  of  the  structure  even  when  excited  at  a resonant 
frequency.  In  cases  where  two  or  more  resonant  frequencies  are  very  close 
together,  it  is  extremely  difficult  to  excite  only  one  mode  shape  of  the  struc- 
ture. The  technique  used  to  calculate  the  structure's  mode  shapes  in  this 
chapter  assumes  that  \ mode  shapes  are  participating  in  the  structure's 
response  over  the  entire  frequency  range  of  interest. 

Finally,  when  attempting  to  excite  only  one  mode  shape  of  a structure 
one  must  excite  the  structure  at  a resonant  frequency.  The  admittances  dealt 
with  in  this  section  are  not  measured  at  resonant  frequencies  of  the  structure 

where  the  dynamic  response  of  the  structure  is  governed  to  a large  degree  by 

1 2 

damping  and  affected  by  shakers  and  suspension  systems.  These  factors 
can  distort  the  relative  magnitudes  and  phases  of  the  measured  components 
of  the  structure's  mode  shapes. 
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CHAPTER  III 


EXPERIMENTAL  GENERATION 
OF  MOMENT -ROTATIONAL  VELOCITY  IMPEDANCES 

When  assembling  structures  into  a total  structure  using  mechanical  im- 
pedance techniques,  moment-rotational  velocity  impedances  are  required  to 
ensure  the  compatibility  of  rotational  displacements  at  the  interfaces  of  the 
constituent  structures.  Usually  it  is  impractical  to  measure  impedances 
directly  from  the  structure  due  to  the  necessary  boundary  conditions  that  all 
DOF  except  one  must  have  zero  velocity  (displacement).  It  follows  that  one 
should  calculate  the  impedances  from  admittances  measured  on  the  structure. 

To  calculate  moment-rotational  velocity  impedances  it  is  necessary  to 
generate  moment -rotational  velocity  admittances.  In  many  situations,  meas- 
ured moment-rotational  velocity  admittances  have  eluded  engineers  due  to  the 
difficulty  in  fabricating  a shaker  that  applies  a pure  moment  to  a structure. 

The  DE  principle  can  be  used  to  circumvent  this  problem.  The  structure 
is  assumed  to  have  X mode  shapes  and  L DOF,  L > X.  The  structure's  DOF 
are  divided  into  two  subsets,  X linear  DOF  in  the  first  subset  and  L-X  rotation 
DOF  in  the  second  subset.  Equation  (25)  can  be  used  to  calculate  moment- 
rotational  velocity  admittances,  [Y,,J,  from  force -velocity  admittances,  [Y^], 
and  force -rotational  velocity  admittances,  [Y^].  Then  one  could  calculate 
the  moment-rotational  velocity  impedances,  [Z^],  fromlY,,]. 

This  method  of  calculating  the  moment-rotational  velocity  impedances, 
although  straightforward,  has  many  shortcomings.  First,  one  must  be  able 
to  measure  accurately  all  the  admittances,  [Y^]  and  [Y^],  over  the  entire 
frequency  range  of  interest.  As  seen  in  the  last  chapter,  admittances  can  be 
distorted  near  resonant  frequencies  by  damping.  Other  errors  can  be  intro- 
duced near  resonant  frequencies  by  shakers  and  suspension  systems.  In 
addition,  the  frequency  range  of  the  measuring  equipment  may  not  cover  the 
entire  frequency  range  of  interest.  Second,  if  the  order  of  [ Y^ ^ greater 
than  X,  the  inverse  of  [ Y^ ^),  [Z^  ,],  does  not  exist  because  there  are  only  X 
linearly  independent  mode  shapes  of  the  structure.  The  third  major  short- 
coming is  the  large  amount  of  data  that  must  be  processed  to  calculate  [Z^l 
at  all  the  necessary  frequencies. 
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It  is  more  efficient  to  solve  for  two  invariant  matrices  of  the  structure, 
the  mass  matrix  and  stiffness  matrix,  and  use  them  to  calculate  an  impedance 
matrix  containing  moment -rotational  velocity  impedances: 

[Z]  = ioj  (M]  + (iw)"1  [K]  . (44) 


Initially,  one  calculates  the  mass  matrix  and  stiffness  matrix  using  the  inverse 

eigenvalue  method  from  the  admittances  measured  at  the  structure's  first  DOF 
1 3 

subset.  Since  the  mass  matrix,  [M^J,  and  the  stiffness  matrix,  [ ^ ] , are 
calculated  from  admittances  measured  only  at  the  subset  of  DOF,  the  mass 
matrix  and  stiffness  matrix  characterize  the  behavior  of  the  structure  only  at 
its  first  DOF  subset. 

The  mode  shapes  at  the  first  DOF  subset,  and  resonant  frequencies 

of  the  structure,  co  , are  calculated  from  the  mass  matrix  and  stiffness  matrix 

1 14 

using  the  iterative  method. 

{0  } = -uf  [Mh)  {<*>*}  + [K^]  {<*>/}  1=1,2,  3...\  (45) 


where 


[<*>/ J IM41]  {<*>/}  = m^ 


{<*>/} 


£ = 1,  2,  3.  . .X  (46) 

£ = 1,  2,  3.  . . X (47) 


= £^  resonant  frequency  in  radians  per  second 
m{  = £ modal  mass 

{<p  = t ^ unnormalized  mode  shape  for  the  first  DOF  subset 


£1^  normalized  mode  shape  for  the  first  DOF  subset. 
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The  normalized  mode  shape  at  the  first  DOF  subset  are  assembled  into  a 
normalized  modal  matrix  of  the  structure: 


(48) 


[<p^/\ jl J , the  normalized  modal  matrix  for  the  first  DOF  subset,  and  [ 1 

and  [ Y^]  measured  in  the  vicinity  of  each  resonant  frequency  are  sufficient 
to  calculate  [02/p],  the  normalized  modal  matrix  for  the  second  DOF  subset 
using  Eq.  (49): 


T 


Equation  (49)  is  derived  from  Eq.  (22)  in  Chap.  I: 


[02]t  = [01)t[y11)-1(y21]t 


(49) 


(22) 


The  normalized  mode  shapes  of  the  structure  are  determined  at  both  DOF 
subsets  and  are  used  to  construct  the  total  normalized  mode  shapes  of  the 
structure: 


! = 1,  2,  3.  . . X 


(50) 


X components,  chosen  from  each  total  normalized  mode  shape,  are  used  to 
construct  each  reduced  mode  shape,  {0^/jjl  }.  The  impedances  to  be  calculated 
from  will  be  used  to  ensure  the  compatibility  of  displacement  at  the 

interfaces  of  the  structure  with  other  structures.  As  a result  the  X compo- 
nents of  must  represent  at  least  all  the  rotational  and  linear  DOFs  on 

the  interfaces  of  the  structure. 

The  reduced  normalized  mode  shapes  of  the  structure,  [0^/p]  and  the 
structure's  resonant  frequencies,  are  used  to  calculate  mass  and  stiffness 
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matrices  representing  the  DOF  at  the  interfaces  of  the  structure. 


2 

1 


0 


co 


2 

3 


= [**£]  • (51) 


0 


co 


2 

X . 


The  relationships  of  the  mass  matrix  and  stiffness  matrix  of  a structure  to  its 
normalized  mode  shapes  and  resonant  frequencies  are  given  in  Eqs.  (52)  and  (53): 


0- 


[K 


3 3 


0. 


0- 


M 


33J 


0- 


= It] 


(52) 

(53) 


T “1 

Premultiplying  both  sides  of  Eqs.  (52)  and  (53)  by  [0,/|ij  and  postmultiplying 
both  sides  of  both  equations  by  [ 0 ^/pc J produces  Eqs.  (54)  and  (55): 


[k33] 


03  T-l 


l^2] 


[0- 


-1 


(M 


T-l 


(54) 

(55) 


Equations  (54)  and  (55)  show  that  the  structure  must  have  at  least  as  many 

linearly  independent  mode  shapes  as  there  are  DOF  at  its  interfaces  or  else 

T™1  “1  T ~ 1 “1 

[0^/V]  and  [(p^/[x]  do  not  exist.  If  [ 0 and[0^/fji]  do  not  exist. 
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Eqs.  (54)  and  (55)  cannot  be  used  to  calculate  the  structure's  mass  matrix  and 
stiffness  matrix.  Using  Eqs.  (54)  and  (55)  one  can  solve  for  the  mass  matrix 
and  stiffness  matrix  of  the  structure  representing  its  DOF  at  its  interfaces. 

Using  the  relationship  in  Eq.  (56)  the  corresponding  impedance  matrix 
can  be  constructed: 

[Z33]  = iu>  [M33]  + (iw)"1  [K33J  . (56) 

Conventional  impedance  techniques  can  now  be  applied  to  assemble  analytically 
this  structure  to  other  structures  ensuring  the  compatibility  of  linear  and  rota- 
tional displacement  at  the  interfaces. 
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CHAPTER  IV 


COUPLING  TWO  BEAM  STRUCTURES 
USING  EXPERIMENTALLY  GENERATED  IMPEDANCES 

The  6 -ft  aluminum  clamped -clamped  beam  structure  dealt  with  in  Chap.  II 
is  cut  in  half.  The  two  halves  of  the  6 -ft  beam  structure  are  reassembled 
analytically  using  impedances  generated  experimentally  as  shown  in  Chap.  III. 
The  analytically  reassembled  beam  structure's  resonant  frequencies  and  mode 
shapes  are  compared  with  the  resonant  frequencies  and  mode  shapes  generated 
for  the  6 —ft  beam  structure  in  Chap.  II.  The  mode  shapes  and  resonant  fre- 
quencies of  a similar  continuum  beam  are  also  shown  for  comparison. 

Both  halves  of  the  6 -ft  beam  structure  are  modeled  as  clamped -free  beam 
structures,  each  with  three  DOF,  L = 3 (Fig.  7). 


18-7-2722 


DOF  DOF 


Fig.  7.  Two  clamped-free  beam  structures. 


Admittances  are  measured  with  equipment  (Fig.  2)  described  in  Chap.  II 
at  the  DOF  of  structure  1 in  Fig.  7.  (Due  to  the  similarity  of  structures  1 and 
2,  admittance  measurements  are  made  only  on  structure  1.  The  impedances 
of  structure  1 are  determined  as  shown  in  Chap.  Ill;  impedances  for  structure 
2 are  constructed  from  impedances  of  structure  1 by  inspection.)  The  actual 
clamped-free  beam  with  shaker/ impedance  head  attached  (Fig.  8)  is  separated 
from  its  support  by  a steel  plate  at  the  base  of  the  clamped-free  beam.  The 
beam  tilts  to  the  right  to  provide  room  for  the  accelerometer  that  measures 
the  response  at  the  rotational  DOF. 
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Fig.  8.  Three-ft,  clamped-free 
beam  structure  with  shaker/ 
impedance  head  attached. 


Fig.  9.  Lever  arm,  accelerometer 
and  shaker/impedance  head  meas- 
uring cross -admittances  for  the 
rotational  DOF  of  structure  1. 


I 


^202-5463 
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The  dynamic  response  at  the  rotational  DOF  of  the  structure  is  measured 


by  an  accelerometer  mounted  on  a lever  arm  to  detect  changes  in  the  slope  of 
the  free  end  of  the  beam  (Fig.  9).  For  small  rotations  the  rotation  of  the  free 
end  of  the  beam  in  radians  is  equal  to  the  slope  at  the  end  of  the  beam: 

& = e 

dx 

(57) 

£ = e 

dx 

(58) 

y _ a 

R “ 0 

(59) 

± = k 

m f 

(60) 

Y/R  = Yq  . 

(61) 

Equation  (61)  shows  that  the  force-linear  velocity  cross -admittance  measured 
at  DOF  3 divided  by  the  length  of  the  lever  arm,  R,  is  equal  to  the  force - 
rotational  velocity  admittance  of  the  free  end  of  the  beam. 

This  relationship,  Eq.  (61),  is  invalid  if  there  is  any  significant  axial  motion 
of  the  beam.  The  shaker/impedance  head  is  attached  only  at  DOFs  1 and  2 
making  the  applied  forces  orthogonal  to  the  axial  modes  of  the  beam.  In  addi- 
tion, the  locations  and  direction  of  the  applied  forces  are  chosen  such  that 
torsional  motion  of  the  beam  is  not  produced. 

The  low  frequency  dynamic  behavior  of  structure  1,  less  than  200  Hz,  is 
dominated  by  two  mode  shapes.  These  mode  shapes  correspond  to  two  resonant 
frequencies  below  200  Hz.  The  third  resonant  frequency  of  the  structure  is 
approximately  500  Hz.  The  modal  participation  coefficient  of  the  third  mode 
is  considered  small  compared  to  the  modal  participation  coefficients  of  the 
first  two  modes  for  frequencies  of  excitation  less  than  200  Hz.  For  low  fre- 
quencies of  excitation  the  structure  is  assumed  to  have  only  two  resonant  fre- 
quencies and  two  mode  shapes,  X = 2,  while  the  structure  has  three  DOF, 

L = 3. 

The  first  DOF  subset  is  DOFs  1 and  2 and  the  second  subset  is  DOF  3. 

The  mutual  and  cross -admittances  measured  at  the  first  DOF  subset  are  used 
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1 5 

to  calculate  the  structure's  mass  matrix  and  stiffness  matrix.  This  mass 
matrix  and  stiffness  matrix  characterize  the  behavior  of  the  structure  only  at 
its  first  DOF  subset: 


[Kll)l  = 


2.96  -1.03 

-1.03  0.43  J 


10  lb/ in. 


(62) 


[Mll)l  = 


' 4.71 
. 0.65 


0.70  ' 
1.70. 


slug 


(63) 


The  structure's  resonant  frequencies,  modal  masses  and  normalized  mode 

shapes  for  the  first  DOF  subset  are  calculated  from  the  mass  and  stiffness 

16 

matrices  using  the  iterative  method.  The  results  are  given  in  Table  3. 


TABLE  3 

RESONANT  FREQUENCIES,  MODAL  MASSES 
AND  NORMALIZED  MODE  SHAPES  FOR  STRUCTURE  1 


The  normalized  mode  shapes  at  the  first  DOF  subset,  [0^/p],  and  the  admit- 
tance matrices  [Y^]  and  [Y^]  measured  at  40  Hz  and  Eq.  (49)  are  used  to 
solve  for  structure  l's  normalized  mode  shapes  at  the  second  DOF  subset. 
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i^ii 
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,-1 
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(64) 


1*2 

T 

6.9 

1 8. 8" 

'-10.86 

4.36' 

'3.25' 

. -2  . 

'-0.585' 

= 

i 

10  i = 

H* 

1 

.-13.4 

16.3. 

. 4.36 

-1.33. 

.6.01 . 

.-2.23 

The  reduced  mode  shapes,  {$^/p}f  are  constructed  using  components  from 
DOFs  2 and  3 and  are  shown  in  Table  4. 

TABLE  4 

REDUCED  MODE  SHAPES  OF  STRUCTURE  1 


18.8 

\tl\ 

16.3 

hi 

1 

-0.585 

P2 

1 

-2.23 

187-2723 


Both  DOF  at  the  interface  of  structure  1 
with  structure  2 are  accounted  for  in 
the  reduced  mode  shape.  As  a result, 
compatibility  of  lateral  and  rotational 
displacement  at  the  interface  is  as- 
sured when  the  two  structures  are 
assembled  using  impedance  techniques. 

The  reduced  mode  shapes,  [ d>^/ p]^  , 
the  structure's  resonant  frequencies 
and  modal  masses  calculated  above  are 

used  to  calculate  the  mass  and  stiffness  matrices  of  structure  1 (Fig.  10). 
mass  and  stiffness  matrices  are  calculated  using  Eqs.  (54)  and  (55)  and  are 
shown  in  Eqs.  (65)  and  (66). 


DOF 


Fig.  10.  Structure  1 using 
DOFs  2 and  3. 
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h 

p 


T-l 


-1 


(54) 

(55) 
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(K33«!  = 


0.466  12.1 

12.1  374 


io- 


(65) 


[m33]1  = 


4.73  44.41 

.44.4  586 


10 


-3 


(66) 


-3 

The  mass  matrix,  Eq.  (66),  is  corrected  for  the  0.27  x 10  slug  mass  of  the 
lever  arm  in  Fig.  9.  Its  rotary  inertia  is  negligible. 

The  mass  and  stiffness  matrices  are  used  to  construct  the  impedance 
matrix  of  the  structure  in  Fig.  10: 

[Z33]  = iw  [M33]  + (iw)"1  [K33]  . (56) 

By  changing  the  sign  of  the  appropriate  elements  in  [ 3 ] ^ and  one  can 

produce  the  mass  matrix,  stiffness  matrix  and  impedance  matrix  of  structure  2 
(Fig.  11). 
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Fig.  11.  Structure  2 modeled 
with  two  DOF. 


The  mass  matrix,  [M~3]_,  and  the  stiffness  matrix,  of  structure  2 

are  given  in  Eqs.  (67)  and  (68): 


[K 


33^2 


0.466 

.-12.1 


-12. r 

374 


(67) 


[M 


33^2 


4.73 


-44.4 


-44.4 

586 


(68) 


30 


(69) 


[Z33)2  iu  [M33J2  + (lU))  [ K3 3 ] 2 

Using  conventional  impedance  techniques  the  structure  in  Fig.  10  is  coupled 
to  the  structure  in  Fig.  11.  The  resultant  structure  is  a 6 -ft  clamped-clamped 
beam  structure.  The  impedance  matrix  of  the  total  structure  [Z]^,  is  deter- 
mined using  Eq.  (70): 

[Z]rp  = t ^ 3 3 1 1 + ^33^2  = [ I ^3  3^1  + l ^3  3^2 

+ (io;)"1|[K33]1  +[K33]2]  (70) 


[ kjt 


[k33j1  +[k*j 


33  J2 


' 0.932 
. 0 


748  J 


(71) 


[M]T  = [M33]1  + [M33]2 


0.946 

0 


0 

117 
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(72) 


Fig.  12.  Total  beam  structure. 

The  total  structure's  resonant  frequencies  and  normalized  mode  shapes 
are  calculated  from  the  mass  matrix,  [M],p,  and  stiffness  matrix,  [K],p.  The 
resonant  frequencies  of  the  6 -ft  clamped  beam  structure  calculated  in  this 
chapter  are  compared  in  Table  5 to  the  resonant  peaks  observed  in  the  admit- 
tances measured  for  Chap.  II,  the  resonant  frequencies  calculated  in  Chap.  II 
and  the  resonant  frequencies  of  a 6 -ft  clamped-clamped  continuum  beam. 
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TABLE  5 


COMPARISON  OF  RESONANT  FREQUENCIES 
FOR  6 -FT  DOUBLE-CLAMPED  BEAMS 


Observed  Resonant 

= 50.0 

Hz 

Continuum 

= 55.3  Hz 

Peaks  (Chap.  II): 

1 

P 2 

= 138.0 

Hz 

Beam: 

1 

"2 

= 152.0  Hz 

Calculated  in 

= 53.2 

Hz 

Calculated  in 

= 50.0  Hz 

Chap.  II: 

1 

V. 

= 139.0 

Hz 

Chap.  IV: 

1 
V ^ 

= 127.0  Hz 

2 

2 

The  second  resonant  frequency  calculated  from  the  experimentally  gener- 
ated impedances  is  relatively  inaccurate  compared  to  the  values  determined 
for  the  second  resonant  frequency  in  Chap.  II.  The  accuracy  of  the  second 
calculated  resonant  frequency  might  be  improved  if  a more  accurate  component 
of  the  second-mode  shape  of  structure  1 at  DOF  3 were  calculated.  The  com- 
ponent of  the  second -mode  shape  of  structure  1 at  DOF  3 is  calculated  using 
Eq.  (49).  The  admittance  matrices  [ ^ ] and  [ ] used  in  Eq.  49  are  measured 

at  40  Hz,  which  is  relatively  close  to  the  first  resonant  frequency  and  relatively 
far  from  the  second  resonant  frequency  of  structure  1.  As  a result  the  partici- 
pation of  the  second  mode  is  far  less  than  the  participation  of  the  first  mode  in 
[ Y^]  and  [Y^^J.  Under  these  conditions  one  would  expect  to  calculate  the 
components  of  the  second -mode  shape  with  less  accuracy  than  the  components 
of  the  first-mode  shape.  It  is  not  possible  to  calculate  the  component  of  the 
second-mode  shape  of  structure  1 at  DOF  3 with  data  measured  closer  to  the 
second  resonant  frequency  due  to  significant  elastic  motion  of  the  lever  arm 
(Fig.  9)  at  frequencies  above  40  Hz. 

If  resonant  frequencies  calculated  for  the  total  structure  are  not  in  the 
frequency  range  where  the  X mode  shapes  dominate  structure  1 or  structure 
2's  response,  the  calculated  resonant  frequencies  and  corresponding  mode 
shapes  are  probably  inaccurate.  The  assumption  at  the  beginning  of  this 
chapter,  that  the  dynamic  behavior  of  structures  1 and  2 is  dominated  by  their 
first  two  modes,  is  not  valid  at  frequencies  above  200  Hz.  If  resonant  fre- 
quencies and  mode  shapes  of  the  total  structure  above  200  Hz  are  desired. 
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then  one  must  characterize  structures  1 and  2 with  additional  mode  shapes 
and  DOF  to  accurately  model  their  dynamic  behavior  above  200  Hz. 

The  mass  matrix,  and  stiffness  matrix,  [K]^,,  of  the  total  structure, 

Eqs.  (71)  and  (72),  are  relatively  accurate  compared  to  the  mass  matrix  and 
stiffness  matrix  of  a 6 -ft  clamped  continuum  beam,  Eqs.  (73)  and  (74).  (The 
methods  used  to  construct  Eqs.  (73)  and  (74)  are  discussed  in  the  Appendix.) 
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The  normalized  mode  shapes  of  the  total  structure  calculated  from  the 
total  structure's  mass  matrix,  and  stiffness  matrix,  [K]^,  are  given  in 

Table  6. 


TABLE  6 

NORMALIZED  MODE  SHAPES  OF  THE  TOTAL  STRUCTURE 
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10.3 

</>2 
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**1  T 
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^2 
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Unfortunately  the  mode  shapes  in  Table  6 include  only  the  DOF  of  the  total 
structure  encompassed  by  the  reduced  mode  shapes  of  the  constituent  substruc- 
tures. These  DOF  are  shown  in  Fig.  12.  The  DOF  at  midspan  of  structure  1 
and  structure  2 (Fig.  7)  are  not  accounted  for  in  the  mode  shapes  in  Table  6. 
These  interior  DOF  are  not  included  in  the  reduced  mode  shapes  of  their  re- 
spective substructures. 

Assume  that  the  total  structure  is  in  steady-state  free  (no  applied  forces) 
motion  proportional  to  only  one  of  its  mode  shapes: 


{x} 


n 


(75) 


The  displacements  at  the  boundaries  of  each  substructure  are  known. 
Focusing  attention  on  one  substructure,  Eq.  (75)  is  rewritten  to  include  only 
the  DOF  at  the  substructure's  boundaries.  These  DOF  are  also  the  DOF  in- 
cluded in  the  reduced  mode  shapes  of  the  substructure,  Eq.  (76): 


n 


(x3}  = c I tr 


(76) 


The  displacements  on  the  boundary  of  the  substructure,  {x^},  can  also  be 
expressed  as  a linear  combination  of  the  substructure’s  mode  shapes: 


(77) 


Since  there  are  X DOF  on  the  boundary  of  the  substructure  where  the 
displacement  is  known  and  there  are  X,  linearly  independent,  reduced -mode 
shapes  for  the  substructure,  there  exists  a unique  solution,  {a},  to  Eq.  (77): 


03.-l 


{x,}  = {a} 


(78) 


Using  Eqs.  (76)  and  (78),  one  can  relate  the  components  of  the  total  structure's 
mode  shape  at  the  boundaries  of  the  substructure  to  the  reduced  mode  shapes 
of  the  substructure: 
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e 


(79) 


All  the  quantities  on  the  left-hand  side  of  Eq.  (79)  are  known,  therefore  (a)  is 
determined. 

Since  the  modal  participation  of  the  reduced -mode  shapes  {a}  of  the  sub- 
structure is  determined  and  the  components  of  the  substructure's  mode  shapes 
at  its  interior  DOF  were  previously  determined,  Table  3,  one  can  solve  for 
the  motion  of  the  substructure  at  its  interior  DOF: 


The  initial  assumption  for  this  derivation  is  that  the  total  structure's  motion 
is  proportional  to  one  of  its  mode  shapes.  Therefore  {x^}  is  proportional  to  the 
total  structure's  mode  shape  at  the  DOF  in  the  interior  of  the  substructure: 


Using  Eqs.  (79),  (80)  and  (81)  the  components  of  the  total  structure's  mode 
at  the  DOF  in  the  interior  of  the  substructure  are  determined: 


The  complete  mode  shapes  for  the  total  structure  shown  in  Fig.  13  are 
calculated  and  given  in  Table  7. 
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Fig.  13.  Total  beam  structure  with  all  DOF. 
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TABLE  7 


COMPLETE  MODE  SHAPES 
FOR  THE  TOTAL  BEAM  STRUCTURE 


For  comparison,  the  normalized  mode  shapes  of  the  6 -ft  clamped-clamped 
beam  structure  determined  in  Chap.  II  are  given  in  Table  8.  Table  9 shows 
the  normalized  mode  shapes  of  a 6 -ft  clamped-clamped  continuum  beam. 


TABLE  8 

NORMALIZED  MODE  SHAPES  FOR  THE  6 -FT  BEAM 
STRUCTURE  DETERMINED  IN  CHAPTER  II 


TABLE  9 

NORMALIZED  MODE  SHAPES  FOR  A 6 -FT 
CLAMPED-CLAMPED  CONTINUUM  BEAM 
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OVERVIEW  AND  CONCLUSIONS 


The  principle  of  dynamic  equipollency  coupled  with  the  inverse  eigenvalue 
method  can  overcome  two  weak  areas  of  mechanical  admittance/ impedance 
testing.  The  principle  makes  possible  the  generation  of  moment-rotational 
velocity  admittances  and  impedances  from  experimental  data  as  well  as  the 
efficient  treatment  of  structures  having  many  DOF. 

When  solving  for  a mode  shape  of  a structure  having  many  DOF  it  is 
necessary  to  take  admittance  measurements  relatively  close  to  the  resonant 
frequency  corresponding  to  the  mode  shape  in  question.  This  ensures  a 
relatively  large  participation  of  the  mode  shape  in  the  measured  admittances. 
It  should  be  kept  in  mind,  however,  that  admittance  measurements  taken  too 
close  to  a resonant  frequency  may  be  distorted  by  nonlinear  damping,  the 
effects  of  the  shaker  and/or  the  effects  of  a suspension  system. 

The  experimental  generation  of  moment-rotational  velocity  admittances 
and  impedances  makes  possible  the  assembling  of  substructures  into  a total 
structure  ensuring  the  compatibility  of  rotational  as  well  as  lateral  displace- 
ments at  the  interfaces.  The  only  significant  restraint  on  using  the  DE  prin- 
ciple to  accomplish  this  is  that  each  substructure  must  have  at  least  as  many 
mode  shapes  (elastic  and  rigid  body)  as  it  has  DOFs  at  its  boundaries. 
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APPENDIX 

CHARACTERISTICS  OF  CONTINUUM  BEAMS 

Throughout  this  report  the  experimentally  determined  characteristics  of 
the  beam  structures  are  compared  to  the  characteristics  of  similar  continuum 
beams.  This  section  deals  with  the  methods  used  to  calculate  the  character- 
istics of  the  continuum  beams. 

17 

The  mode  shapes  of  the  continuum  beams  are  lifted  directly  from  tables 
and  are  normalized  in  a manner  similar  to  the  experimentally  determined  mode 

shapes  with  which  they  are  compared.  The  resonant  frequencies  of  the  con- 

1 8 

tinuum  beams  are  calculated  using  Eq.  (A-l): 

u>n  = /3n2  4Wp  (A-l) 

where 

= nth  resonant  frequency  expressed  in  radians  per  second 

E = Young's  Modulus  for  aluminum  = 9.67  x 10^  psi 

I = second  moment  of  area  of  the  cross  section  of  the  beam 
= 0.25  in."* 

-4 

p = mass  per  unit  length  of  the  beams  = 3.74  x 10  slug/in. 

/3  = parameter  based  on  length  of  the  beam  and  displacement 

boundary  conditions,  i.e.,  clamped-clamped,  clamped- 
19 

free,  etc. 

The  stiffness  matrices  of  the  continuum  beams  are  calculated  using 
Eq.  (A-2): 

[K]  =[£]T‘1  (A-2) 


38 


where 


= diagonal  matrix  containing  the  resonant  frequencies 
corresponding  to  the  mode  shapes  in  [0/ p] 


[ 0/ pi]  = normalized  modal  matrix  of  the  continuum  beam  and 


contains  a finite  number  of  mode  shapes. 


This  represents  the  assumption  that  the  low  frequency  dynamic  behavior  of 
the  continuum  beam  can  be  approximated  by  its  first  two  mode  shapes.  The 
mode  shapes  in  [0/p]  correspond  to  the  A mode  shapes  of  the  beam  structure 
being  compared  to  the  continuum  beam. 

Each  mode  shape  in  [0/p]  is  normalized  by  the  square  root  of  its  respective 


modal  mass.  The  modal  mass  of  the  continuum  beams  are  calculated  using 
Eq.  ( A — 3 ) : 


pi  = m. 


n 


(A-3) 


where 


m = modal  mass 


n 


( = length  of  the  beam 


p = mass  per  unit  length. 

? 0 

The  tabulated  mode  shapes  for  the  continuum  beams  are  normalized 
21 

using  Eq.  (A-4): 


( A-4) 


The  n^1  modal  mass  of  a continuum  beam  is  defined  by  Eq.  (A-5): 


2 

0 p (x)  dx  = m 
n f ' ' n 


(A-5) 
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Since  p(x)  is  a constant  in  this  case,  Eq.  (A-5)  is  re-expressed  as  Eq.  (A-6): 


C*  2 

P \ (Pn  dx  = m . (A-6) 

Substituting  Eq.  (A-4)  into  (A-6)  yields  Eq.  (A-3). 

The  normalized  modal  shapes  of  the  continuum  beam  are  also  used  to 
calculate  the  mass  matrix  of  the  continuum  beam,  Eq.  (A-7): 

[M]  = [^-]T_1  [^-f1  • (A-7) 

r-  r 

As  before,  the  normalized  mode  shapes  in  [#/(i]  correspond  to  the  X mode 
shapes  of  the  beam  structure  being  compared  to  the  continuum  beam. 
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